Invariants from classical field theory 



Rafael Diaz and Lorenzo Leal 



Abstract 



We introduce a method that generates invariant functions from perturbative classical 
field theories depending on external parameters. Applying our methods to several field 
theories such as abelian BF, Chern-Simons and 2-dimensional Yang-Mills theory, we obtain, 
respectively, the linking number for embedded submanifolds in compact varieties, the Gauss' 
and the second Milnor's invariant for links in S^, and invariants under area-preserving 
diffeomorphisms for configurations of immersed planar curves. 

1 Introduction 

Suppose we have a physical system determined by an action functional S depending on fields 
as usual, but also on external parameters p belonging to some space P. Fixing the external 
parameters we can study our field theory from a classical or quantum point of view. A funda- 
mental problem in physics is to understand how the system changes as p varies. In this work 
we focus our attention on a particular issue arising within this setting. Suppose that fields 
and external parameters are acted upon by a Lie group G, and that S is invariant under the 
simultaneous action of G on fields and parameters. How is the group of symmetries G reflected 
on the states of the system as p varies? The answer to this problem is very different depending 
on whether we look at it from a quantum or a classical point of view. The quantum situation 
requires techniques related to anomalies and will not be discussed in this work. 



We look at the classical situation from two points of view: we consider exact as well as 
perturbative solutions of the equations of motion. The exact results admit a simple and clean 
description given in Section [2j Simplifying slightly, the conclusion is that associated with such 
a classical system there is a function Sos '■ P — > 1^ whose value at p is obtained by evaluating 
S on-shell, i.e. Sosip) = S{ip{p),p) where ip{p) represents a solution of the equations of motion 
of the action S{ ,p). If the correspondence p — > ip{p) can be constructed in a G-equivariant 
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fashion, then the function Sqs ■ P — > IK is invariant under the action of G. We apply this tech- 
nique to abelian BF theory on compact oriented manifolds, with the external parameters being 
pairs of non-intersecting embedded submanifolds of the appropriated dimensions. The on-shell 
action is the linking number of the embedded submanifolds, an invariant under diffeomorphism 
of the ambient manifold connected to the identity. We write explicitly using local coordinates 
the on-shell action for this example. 

It is seldom possible to find explicitly the solutions of the equations of motion for most inter- 
esting Lagrangian systems. In this case an alternative route is to look for perturbative solutions 
of the equation of motion. Section |4] contains the main result of this work: the proof that a 
hierarchy of invariant functions on parameter space can be obtained from the computation of 
the perturbative on-shell action Sqs of a functional action S invariant under the simultaneous 
action of a Lie group on fields and external parameters. In order to show this result we de- 
velop in Section [3] a fairly explicit model for the study of perturbative solutions of systems of 
equations (algebraic, differential, integral, etc.) We show that, in the non-degenerate case, the 
perturbative solution of the system is unique. In the degenerated case, even though uniqueness 
is lost, our methods guarantee that the space of perturbative solutions is non-empty, and pro- 
vide an explicit solution under weak assumptions. We apply this technique to solve a couple of 
general equations arising in the context of Hodge algebras using techniques closely related to 
homological perturbation theory [18j. In Section [5] we consider an example of interest in low 
dimensional topology, namely, we applied our methodology to Chern-Simons-Wong action and 
show that it yields link invariants, in particular, we obtain the Gauss' and the second Milnor's 
invariants. In Section [6] we discuss how invariants under area preserving diffeomorsphisms of 
'S? can be obtained given a generic finite family of immersed curves in the plane, applying our 
methods to Yang-Mills- Wong action. The Wong term that we add both in the Chern-Simons 
and Yang-Mills theories may be thought, physically, as the action for conservation of cromo- 
electric charge, and mathematically, as the action functional for parallel transport in a fiber 
bundle. In Section [7] we present a brief discussion of open problems and future lines of research. 



2 Exact results 

Fix spaces F and P, thought as the space of fields and parameters, respectively. Assume that a 
Lie group thought as the gauge group, acts on F. Fix another Lie group G which acts on F 
and P, together with a map k : G — > K. Let G act on F x P via the diagonal action. Suppose 
we have map S : F x P — > M, thought as the action of a classical field theory, satisfying 
S{gip,gp) = S{ip,p) and S{kip,p) = S{ip,p) for {ip,p) ^ F x P, g G G and k £ K. In addition, 
assume we have map a : P — > F such that a{gp) = k{g)ga{p) for g (z G and p £ P. Lemma [T] 
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below explains how one can get a G-invariant function on P from this data. 

Lemma 1. The map Sa ■ P — > given by Sa{p) = S{a{p),p) is G-invariant. 

Indeed for g G and p € P we have that: 

Sa{gp) = S{a{gp),gp) = S{k{g)ga{p), gp) = S{ga{p),p) = S{a{p),p) = Sa{p)- 

A fundamental question thus arises in this context: how can one obtain such a map a? 
We are going to show via examples that it is often possible to find a map a with the required 
properties by solving the equation of motion, i.e. finding for each p € -P a solution ip{p) of the 
equation 

Thus we are going to show that the so-called on-shell action Sos{p) = S{ip{p),p) is a G-invariant 
function on parameter space. 

Let us first consider abelian BF gauge theory generalizing a construction of [24j . Let M be 
a compact oriented manifold of dimension n and fix 1 < p < n. The space of fields is 

where Q^(M) denotes the space of differential i-forms on M. Let BE{M,i) be the space of 
bounding embedded i-dimensional submanifolds of M, i.e. 

BE{M,i) = {7 I 7 : S — > M bounding embedding, S a compact oriented i-manifold }/ ~ . 

An embedding 7 : S — > M is bounding if there exists embedding 5 : A — > M, where A is an 
oriented manifold with boundaries such that d{A) = S and 6\-e = 7. Embeddings 71 : Si — > M 
and 72 : S2 — ^ M are ~ equivalent if there exists an orientation preserving diffeomorphism 
(p : Si — > T,2 such that 71 = o 72. The space of parameters is BE{M,p) x BE{M, n — p — 1) 
and the action functional 

S : (17P(M) eJl"-P-^(M)) X BE{M,p) x BE{M,n-p-l) — >R 

is a BF theory couple to external parameters given by 

-5(^1,^2,71,72)=/ AiAdA2+ 1IA1+ 72*^2, 

for (^1,^2,71,72) e (l^P(M) efi"-P-i(M)) X BE{M,p) x BE{M,n-p-l). The action S is 
invariant under gauge transformations Ai ^ Ai + dfi, A2 A2 + d/2, where /i G QP~'^{M) 
and /2 G il"'~^~^(M). Let A{M) be the infinite dimensional Lie group of automorphisms of M 
connected to the identity; A{M) acts on forms and embedded submanifolds by pull back and 
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push forward, respectively. S is manifestly A(M)-invariant since it is metric- independent. 

Recall [IT] that the Poincare dual form ^(7) G Q'^~p{M) of an embedding 7 : S — > M is 
uniquely determined, modulo de addition of an exact form, by demanding that 

[ YA= [ P{j)AA 
Jm 

for A € 0,P{M). Poincare dual forms have the following properties: if 7 : S — > M is the 
boundary of (5 : A — > M, then d{P{A)) = P{T,); if (j) : M — > M is a diffeomorphism, 
then P{cl)^^ ° 7) = 4'*P{l)- Poincare dual forms, among many other things, are useful to 
compute the linking number //c(7i, 71) of embedded bounding submanifolds 71 : Si — > M and 
72 : T,2 — > M of dimension p and n — p — 1, respectively, as follows: 

^A:(7i,72)= / P(Si)AP(A2) 
Jm 

where 82 ■ A2 — > M is such that d{A2) = S2 and d{62) = 72- Using Poincare dual forms the 
action S may be written as: 

-5(^1,^2,71,72) = / AiAdA2+ [ P(Si)A^i+ / P(S2)A^2. 
Jm Jm Jm 

Varying S with respect to Ai and A2 we obtain the equations of motion 

dAi = (-1)PP(S2) and dA2 = (-lf("-P)+ip(Si). 
Thus the on-shell action —£'05(^1,^2,71,72) is given by 

/ P{^2)/\A2= [ 72*^2= / 72*^2= / 72('^^2)= / P(Si)AP(A2). 

Jm Jt.2 Jd{A2) JA2 Jm 

We have shown that the on-shell action is given by 

-5os(7l,72)= / P(Si)AP(A2). 

Jm 

Prom this expression it is clear that the on-shell action Sqs is an yl(M)-invariant function on 
BE{M,p) X BE{M,n — p — l), indeed 505(71,72) computes the linking number of the embedded 
submanifolds 71, 72- Let us consider the case where M = and use coordinates (xi, X2, ■■■,Xn) 
to write the solution of the equations of motion and the on-shell action. The Poincare dual 
form P(7) = P{'y)ij,dXfj_ of an embedded p-manifold 7 : T, — > is given by 

P{l),{x) = e,,. ^ ^*{dxnS^{x - 7(a)), 
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where a G S, 6"'{x — 7(a)) is the n-dimensional Dh'ac's delta function centered at 7(a), and for 
/i = {fii, ...,/Up) we set dx^ = dx^^ A ... A dx^". The solution = A^^^dx'^ of the equations of 
motion is given by 

(x - vY- 



A^,^(x) = / dy^e^c. 7 ^P(72).(y), 

where |x| = y^x^ + X2 + ... + x^ and e^cv = £iJ,i^i2...iipcui..M„^p^i is the completely antisymmetric 
symbol in n-dimensions. Using the expression above for both Ai and A2, we obtain explicitly 
the on-shell action Sqs, i-e., the linking number of the embeddings 71 and 72: 

Q ^ \ f f *M (71(a) -71 

For n = 3, p = 1, 7i and 72 are actually closed curves in and the on-shell action 
"5*05(71572) is the Gauss' linking number. We may also be consider the case n = 1, p = 0, 
indeed let us consider a slightly generalized action. The space of fields C°°(M)" consists of n- 
tuples (/i, fn) of piecewise smooth functions on the real line M. The manifold of parameters 
Cn(K) = {{xi, x„) E E" I Xi 7^ Xj for i 7^ j} is the space of configurations of n distinguish- 
able points on the real line. Configuration space and its compactification are studied in [UIH]. 
The action 5 : C~(M)''" x Cn(M) — > R given by 

5(/i,..,/n;xi,..,x„) = ^ / fi{x)fj{x)dx -^^fiixi) 

is invariant under under the natural action of the group A{R) of orientation preserving diffeo- 
morphisms of M. The equations of motion are — /j+Z^i<j fj~^Xi = 0, where 5x is the delta 
function concentrated in x and 1 < i < n. Integrating we get — Ylj<i fj + X^«<j fj ~ ^^t-- where 
Ox denotes the Heaviside theta function with jump at x. Thus the equation of motion is Af = 9, 
where A is the n x n matrix given Aij = sg{j — i), f = (/i, ■■■,fn) and & = (^xn •••)^x„)- One 
can check that / = B6 where Bij = (— l)'*^-'' and that the on-shell action Sqs ■ Cn{R) — > R 
is given by 

Sos{xi,...,Xn)= J2 (-1)"""'"'"''^=.'.(^A.). 
i<j,k,s 

Notice that Sosixi, ...,x„) is indeed an j4(M)-invariant function on configuration space as it 
should according to Lemma [TJ 



3 Perturbative solutions 

Suppose one is interested in finding solutions of an equation of the form 0{ip) = ip, where V is 
a vector space, O : V — > 1^ is a non-necessarily linear map, ip is an element of V, and (p is the 
unknown. We are actually going to work perturbatively, so we may as well start with a map 
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O : V — > thus O = O"^" where 0„ : V — F is a non-necessarily linear map. 

Assume that each 0„ admits a globally converging Taylor expansion 

oo 

On = '^On,k{V, ...,9'), 
k=l 

where On,k '■ V^'' — > is a multilinear operator and On,i = for n > 1. Finding solutions 
of the equation Yl'^=o On^^ = V is a notoriously difficult problem, and no general answer 
should be expected. Remarkably, it can be treated perturbatively as follows: performing the 
substitutions 99 — > A(/? and V' ^ip, the equation X^^o On^^ = ip becomes 

^ 0„,fc((^,...,<^)A"+'=-i=V. 

n>0,A:>l 

Making the expansion ip = Y2i=ofi'^^ transforms the equation above, a system of non-linear 
equations, into a system of infinitely many linear equations. Indeed, taking into account the 
powers of A we get a zero order equation: 

Oo,i(¥'o)=^, (1) 
and for n > 1 we get higher order equations: 

Om,ki^il, ■■■,^ik) = 

m,k,ii,...,ik 

where the sum runs over non-negative integers m, k, ii , ij, such that n = m + X^g^^ ig + k — l. 
The equation above may be written in the suggestive form 

Oo,l(v3n) = - X] 0.m,k{^h,-,^i^), (2) 
m,k,ii ,...,ik 

where m > 0, > 2 and n = m-\- "^2^=1 is + k — 1. Thus necessarily integers ii, ...,ik are strictly 
less than n, and if Oo,i is invertible then ([2]) uniquely determines in terms of ipi with i < n. 

In order to find explicitly we need several combinatorial notions A directed graph 
is a triple {V,E,{s,t)) where V and E are finite sets - the set of vertices and edges - and 
(s, t) : E — > V xV is a map. A path 7 in a graph is a sequence of edges ei, 62, such that 
t{ei) = s(ei+i) for 1 < i < A; — 1. We say that 7 is path from s(ei) to t(efc). A rooted tree T 
is a directed graph with a distinguished vertex r, called the root, such that for each vertex v 
of T there is a unique path in T from v to r. The valence of a vertex v is val{v) = \star{v)\, 
where star{v) = {e (z E \ t{e) = v}. Vertex t> is a leave if val{v) = 0. A vertex that is not 
a leave is called internal. The set of internal vertices is denoted by Vj and the set of leaves is 
denoted by Vl. The root of a tree is an internal vertex, except in the case of the tree • whose 
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unique vertex is the root. A tree together with a hnear order on star{v), for each internal 
vertex v is called planar. A labelled planar rooted tree is a pair (T, /) where T is a planar 
rooted tree and / : Vj{T) — > N is a map, called the labelling of T. Consider the category whose 
objects are labelled planar rooted trees. A morphisms is a pair {f,g) where / : Vti — > and 
g : — > are maps such that (52,^2) ° 9 = {f, f) ° isi,ti); moreover we demand the pair 
{f,g) preserves both the label and the linear ordering on star{v) for each internal vertex v. 

Definition 2. For n > 1 let be the set of isomorphism classes of labelled planar rooted 
trees T such that: val{v) > 2 for v & and X^^gV/ (''^"K^) + K^)) = ^ + 

A labelled planar rooted tree T is uniquely constructed by joining planar subtrees Ti, Tk, 
k >2, to the root r labelled by I, see FigurelH If T is so constructed we write T = (Ti, ...,Tk)i. 
The following set theoretical identities hold: 

Vl{T) = u';^MTs) and ^/(T) - {r} = uJ=iF,(r,). 




Lemma 3. Let Tg belong to Tj^ for 1 < s < k. Then T = (Ti, ...,Tk)i belongs to T„ for 
n = Yl''s=i^s + k + l -1. 

Proof. The set theoretical identities above imply that 

{val{v) + l{v)) = ^ {val{v) + l{v)) + valirr) + I = 
veViiT) s,veViiTs) 

^is + Yl + k + l = Y.is + k + l-l + \Vi{T)\. 

s=l s=l s=l 

□ 

Next definition assumes that the operator Oo,i is invertible. 
Definition 4. For T € T„ let Ot '■ V'^ l^^-C^)! — > V be recursively given by 
0. = 0,l and O^T^,...,n), = -Ooi{OiMOT„:;OT,)). 
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Proposition 5. The perturbative solution ip = Y1^=q V^n-^" of equations ([H) and ([2]) is given 
by 

^n=Yl Ot(V',-,V'). (3) 



tgt„ 



Proof. Let ipn be given by ([SD then Oo,i((/?n) = X]tgt„ ^o,i(Ot(^, V'))- By Definition 
Lemma [3l induction and writing T = (Ti, ...,Tk)i, the previous sum is equal to 



^ -0„,fc(OTi(V',---,V'),-",C>Tfe(V',---,^)) = X] -0„,fc(v3i,, ...,(/?iJ. 

n>0,A:>2,TiGTi^,...,TfeeTi^ n>0,fc>2,ji,...,ife 

Thus ifn satisfies the required recursion. □ 

Consider the polynomial equation a„j;" + ... + a2X^ + aix = y with ai 7^ 0. Instead of 
looking for an exact expression for a; as a function of y we look for a perturbative solution, 
i.e. a solution x = X^J^Lo^"-^" °^ ^^e equation OnX'^X^'^ + ... + a2x'^X + aix = y. Let T° be 
the set of isomorphism classes of rooted planar trees with as the label of all internal vertices. 
Proposition [5] implies that 

where 

v(LV^(T) 

Corollary 6. 2; = 'Yl'n=o I'^nl-^" ^^e formal solution of X]^2 — x = —\. 

Similarly one can check that: 

Corollary 7. x = Yl^=o [TnlA" is the formal solution of Ylm k>2^"^"~^''~^ — x = —1. 

We say that an operator O : V — > V has a right inverse if there exists an operator 
P : V — > V such that 0{P{ip)) = for 9? G 0(y\ The proof of Proposition [5] yields the 
following result. 

Proposition 8. Let us assume that "0 € Oi(M), Oo,i posses a right inverse P, and that 
En>o,fc>2,ii,...,ifc ^^^^(v^j^, G Oi(y), where n = m + Es=i^* + A; - 1. A solution of ^ 

and Q is given by = ^^TgT,, OtIV', "0), where ©(j-^^^^j^^), = -P(0;,fe(OTi , OtJ), 
0, = P. 

We refer to the conditions of Proposition [8] as the consistency conditions. To illustrate how 
Proposition [8] works we solve perturbatively two general equations arising in Hodge algebras 
[m [22I [M] . The methods we use to solve these equations resemble the techniques of homological 
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perturbation theory [18]. Let {A, d,<, >) be a Hodge algebra, i.e. {A, d) is a differential graded 
algebra, < , >: Ai^A — > M is a graded symmetric non-degenerated bilinear form, and A admits 
a Hodge decomposition. The adjoint d* of d is such that < da, b >=< a, d*b > for a,b & A. Let 
A = dd* + d*d be the Laplace-Beltrami operator. The subspace 7i A of harmonic elements 
is TC = Ker{A). There is an orthogonal decomposition A = Im{d) ^ Ini{d*) ^Tl, and by 
Hodge theory W is canonically isomorphic to H^{A) = l^^jt-}^ ■ Also there exists an operator 
Q : A — > A such that / = AQ + tt-j-i, where / : A — > A is the identity map and n-n : A — > TC 
is the orthogonal projection onto H. Moreover setting G = d*Q we get that / = Gd + dG + ny,. 
We look for a perturbative solution of the equation 

A(a)+ Yl 0„,fc(a,...,a)A"+^-i =^ (4) 

n>0,k>2 

where On,k ■ A'^^ — > j4 is a linear operator of degree 1 — k, On,i = for n > 1 and b G A^. If 
Ti.^ = 0, then AQ = I on A^ and the consistency conditions of Proposition [8] hold. 

Proposition 9. If TC^ = 0, then a perturbative solution a = X^n^o & given by 

On = Etgt„ Ot(^.-,&) where 0(^Ti,...,n)i = -QiOi,k{OT^, ...,OTk)) and O, = Q. 

Next we look for a perturbative solution of an equation of the form 

da+ Yl 0„,fc(a,...,a)A"+'=-i =6, (5) 

n>0,fc>2 

with a e A^, b G A^, On,k ■ A®^ A operators of degree 2 - k, On,i = for n > 1, and 
d{b) = 0. Assume that the operators On,k satisfy the generalized Leibnitz rule 

k 

dOn,k{ai, Ofc) = ^(-l)^+-+^0„,fc(ai, d{ai), ...,ak), 

i=l 

for homogeneous elements ai, S ^4, where a denotes de degree of an homogeneous element 
a £ A. Moreover, assume that the operator On.k satisfy, for t > 1, the quadratic relations for 
fixed 71, m, t: 

Y Yl i-l)'^^-^~On,kiai,...,Om,iiai,-,ai+i-i),...,at) = 0. (6) 

k+l=t+l l<i<t-l+l 

If H"^ ~ if^(y4) = 0, then nn = and I = Gd + dG on and thus G is a right inverse of d. 

Proposition 10. If = 0, then a perturbative solution of (0) is given by a = X^^q ^nX^ 
where a„ = X]tgt„ ^'t (&,•••, &) and O^T^^^^^^Tk)i = -G{Oi^k{OTi, -^Ot^,)) and O. = G. 

Proof. We have to show that Etgt„ dOrib, b) = 0, where 0(^Ti,...,n)i = -0«,fc(OTi , OtJ- 
Since 

k 

d{-Oi,k{OT„.:,OT,))=Y^On,k{OT„:;dGdT,,...,OT,) = 

i=l 
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k 
i=l 

Thus ^j^grp^ dOr (&,•••, &) equals the sum of two terms; the first one vanishes by ([6]) and the 
second one by induction. □ 

Notice the similarity between the conditions of Proposition [10] and the axioms defining ^oo~ 
algebras [iHl [30l [31] , especially when the operators On^k vanish for n> 1. Indeed our conditions 
involve a countable family of operators O^^k satisfying a countable number of quadratic equa- 
tions. It would be interesting to investigate the operadic and geometric interpretation of the 
conditions of Proposition [TUJ 



4 Perturbative on-shell action 

A major difficulty in the process of obtaining invariant functions by evaluating the on-shell 
action of classical field theories is that one can seldom find explicitly the solutions of the equa- 
tions of motion. We show in this section that one can get around this problem if we evaluate 
instead the perturbative on-shell action. An interesting feature of the perturbative approach 
is that one gets automatically a hierarchy of invariants indexed by the natural numbers. The 
zero level is obtained by linearization of the equations of motion. Higher order invariants are 
obtained applying a sophisticated recursive procedure, where each step consists in solving the 
linear equations of motion with varying non-homogeneous term. 

We are ready to discuss the main result of this work. We are going to show that under 
suitable conditions, made precise below, if we are given an action 

S :F xP — > m[[X]] 

then there are infinitely many G-invariant functions S(^n) '■ P — ^ with n > that are 
constructed by evaluating the on-shell action Sqs perturbatively. Let us then proceed to state 
the conditions necessary for this result. First, we assume that we have a group G which acts 
via a diagonal action on F x P. Second, we assume that the space F of fields is provided with 
a non-degenerated G- invariant symmetric bilinear form < , >: F ^ F — > M. Moreover, we 
assume that linear operators on F can be written in the form < x, > for some x ^ P- Expand 

5 in powers of A as 

oo 

5(v?,p) = ^5„((^,p)A" 

n=0 

and consider the further expansions 

So{f,P) = 2_^ — ■ 7 and for n>l set Sn{(p,p) = — ■ 7 > 
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where Qn,k ■ F®'' x P — > R. Notice that for n > 1 the maps Qn,k ^-re defined for /c > 3, 
this assumption fits nicely with the results of the previous section. Our third assumption is 
that Qn.kio^i 5V) 9P) = Qn,ki^^ ■■■tV,p) for each g G G. By the previous assumptions we can 
write Qn,k{f,---,'P,'4',p) =< On,k-i{'P, ■■■,f,p),'4^ > where the maps On,k ■ F'^'' x P — > F are 
such that 

On,ki9^,-n9^,ap) = gOnM'P^-^'P^p) 

for g & G. Set also Qo,i{ip,p) = — < jip),^^ > where j{gp) = gj{p) for g G G. The Euler- 
Lagrange equations are determined by the identity 

n>0,fc>3 

By the previous assumptions the critical points of S are the solutions of the equation 

Oo,l{^,P)+ ^ On,k{^,:.,(p,p)X'^ = j{p). 
n>0,k>2 

Making ip — > Xip, j — > Aj, the critical points of S are determined by the equation 

Oo,i((p,p)+ 0„,,(<^,...,(/.,p)A"+'^-i=j(p). (7) 

n>0,fc>2 

Proposition 11. If Oi( ,p) is invertible for each p ^ P, then the perturbative solution 
^(p) = E^o'^n(p)A" of dH) is such that ip{gp) = g(p{p) for g £ G. 

Proof. We show that (pn{gp) = gfnip) for g (z G. From Proposition [5] and Lemma [12] we get 
fn{gp) = ^ OT{j{gp),...,j{gp),gp) = g ^ OriJip), ...,j{p),p) = g^n{p)- 

TgT„ TeT„ 



□ 



Lemma 12. OT{ga, gp) = gOxia, /3,p) for g G G and a, (3 G F. 
Proof. Assume that T = (Ti, ...,Tk)i, then 

0(Tu-,Tk)ii9a,...,gP,gp) = -O^liOi^kiOTAga, ....gn), ....Or^igr, ...,gl3),gp),gp), 
which by induction is equal to 

-gOo,l{Oi,k{OTi{a, ...,k), ...,Ot^:{t, ...,P),p),p) = 5-0(^1, „,, t^), ■■■,P,p)- 

Similarly one can prove the following result. 



□ 
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Proposition 13. Assume that j{p) E Oi{V,p), Oo,i( ,p) has a right inverse P( ,p) and 
J2n,k>2,ii,...,ik On,k{^ii^--n^ik^P) ^ Oo,i(^'P) where n = m + Yl''s=i'^s + k - I. The solution 
"Pnip) = Etgt^ OtO', of dZj) satisfies = ipn{p) for € G. 

According to Propositions [T3l if Oo,i( iJ*) has a right inverse, then ip{p) = Yl'^=o^n{p)^^ 
given by ([3]) is a perturbative solution of ([7]). Plugging this solution in 5 we obtain that the 
perturbative on-shell action Sqs ■ P — > ^[[^^l] which is given by 

oo 

n=0 

We proceed to show that the functions Si^n) '■ P — ^ 1^ ^ire G-invariant. For n > let R„ 
be the set of isomorphisms classes of labelled planar rooted trees T that can be written as 
T = (Ti, ...,Tk)i, where Tg G Tj^, Y^gis + I = n, for 1 < s < k and A; > 1 if / = 0, and /c > 2 if 
/ > 1. 

Definition 14. For T G R„ let Qt : F^I^^C^)! — ^ R be given by 

QiTi,...,Tk)i = Ql,k{OTi,---,OT^). 

The proof of the following result is similar to that of Proposition [5j 
Proposition 15. S'(„)(p) = EtgRu QtU(.p)^ ■■■,j{p),P) for n > 0. 

We are finally ready to state and prove the main result of this paper. 
Theorem 16. S'(„) : P — > M is a G-invariant function for n > 0. 
Proof. If p £ P and g £ G then 

S{n){gp) = X] QTU{9P),-,j{9P),9P) = J2 QT{j{p)^--^3{p),P) = S(n){p). 
TgR„ TeR„ 

□ 

5 Chern- Simons- Wong theory and link invariants 

The relation between Chern-Simons theory and link invariants was first study in [32] and is by 
now a solid theory [151 ISOl [271 [28], studied from a variety of points of view. A common feature 
of these approaches is that they work at the quantum level. It was proposed in [22] that it 
is possible to construct link invariants from perturbative classical non-abelian Chern-Simons 
action with an extra term due to Wong [H [33]. Our desire to understand the mathematical 
foundations underlying the methodology of [23j was the primary motivation for this work. The 
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results of this section illustrate the full power of Theorem [TBI which yields a hierarchy of invari- 
ant functions starting from functional actions depending equivariantly on external parameters. 
Let be the unit 3-sphere and g be the Lie algebra of a compact semi-simple Lie group G. 
Fix a symmetric non-degenerated bilinear form Tr on g invariant under the adjoint action. The 
space of fields 

{n\S^)^Q) X M{s\Gr 

consists of tuples {a, gi, gn) where a G 17^(5"^) (8)0 is a g-valued 1-form on S^, and gi : — > G 
is a G-valued map on the unit circle. The space of parameters En{S^, S^) consists of n-tuples 
(71, ...,7n) such that 7^ : — > is a embedded closed curve in S"^, and the images of the 7^ 
are mutually disjoint. Physically, a represents the gauge potential and gi the chromo-electric 
charge of a point-particle undergoing non-abelian interactions. The trajectories in of these 
particles are elements of the parameter space, and we will show that the linking of these parti- 
cles is tested in the process of computing the perturbative on-shell action. 

Let A{S^) be the group of automorphisms of connected to the identity. A{S^) acts 
by pullback on n^{S^) (g) g, trivially on M(S'\G)", and by push-forward on En{S^,S^). To 
construct the action functional we introduce some notation. The pullback of a to via 
7j : — > is denoted by ai{t), where t is the standard coordinate on S^. Fix elements Q G g 
and for gi G M{S^,G) for 1 < i < n, let the chromo-electric charge be Ci{t) = gi{t)cig^^ {t) , 
for t G S*^. The covariant derivative of gi along the z-th particle is Dtgi = dtgi + Xai{t)gi. The 
action functional is 

5'(a,ffi, ...,5n,7i, •••>7n) = / Tr{a Ada+'^Xa^) 5*"*(a,5ri, c?„,7i, ...,7„), 

where 

n ,. 

5^"*(a,ffi,...,9n,7i,-,7n) = ^ / dtTr{hgr\t)Dtgi{t)) 

i=l -^T' 

corresponds to the interaction of n classical Wong particles carrying non-abelian charge [H I33j . 
Chern-Simons action is invariant under the group M(5'^, G) of gauge transformations connected 
to the identity. The action of u G M(S'^ G) on a G Q}{S^) (g) is given by o" = u ^au + u ^dt. 
The action S"*"* is gauge invariant if we set c" = Cj and g^ = u~^gi. Non-abelian charges Ci{t) 
transform in the adjoint representation Ci{t)^ = u{t)~^ Ciu{t) . With these conventions S is an 
j4(S''^)-invariant function. 

The variation of S with respect to a yields the equation 

1 " 

Fa = -Y,P{luCi{t)), 

i=l 
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where Fa = da + §[a, a] is the curvature of a and P{'yi,Ci{t)) is a Poincare dual form defined 



2 

via the identity 



[ Tr{P{-fi,Ci{t))Abi)= [ Tr{ci{t)bi{t))dt. 



The variation of S with respect to gi yields the equation DfCi = Ci + A[aj, q] = of conservation 
of non-abclian charges. Thus c,;(t) = Ui{t) c-i u^^{t), where Ui{t) = Pexp ( — A /jj a„;(,s) ds ) is the 
path ordered exponential of the gauge potential a along the curve 7i. According to our general 
theory the on-shell action may be expanded as 

oo 

Sos{H,-,ln) = ^ 'S'(^)(7l,->7n)A" 
m=0 

where each S'(^) should be an A(S'^)-invariant function of the link (71,..., 7„) G En{S^,S^). 
Using gi{t) = Ui{t)gi{0) one can check that 5'*"*(7i, 7„) = 0. Let ^: g — > End{Q) 
be the adjoint representation of q given by x{y) = [x, y] for x,y G g. From the equation 
DtCi = Ci + Xai{ci{t)) = we see that Ci{t) = Pexp{—X jQaidt)ci, and thus the equation 
= ^TH=iP{li,Ci{t)) becomes 



^ 1 " ^ 00 n 

da = --[a,a] + -Y^ P{li)ci + 3 P{luCm,i{t))ci\^ , 

i=l m=l i=l 

where 

» m 

Cm,i{t)= A4j(^)' ^0,m = {i^l^X2,-,Xm) \ < Xj < t and Xj < Xk H j < k}, 

and the map eij : A^^ — ^ is given by eij{xi,X2, ■■■,Xm) = li{xj). 

We look for a perturbative solution a = X^pO(p)A^ of the equation of motion. The cor- 
responding recursive system of linear equations is given by da^o) = 5 Z^"=i -P(7i)ci, and for 
p>l 

1 ^ P n 



si+S2=p— 1 ■m=l i=l 



where 



siH |-Sm=p-m j=l 

Similarly the perturbative on-shell action is Sos = J2m=o '^{m) where for m > we have 
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Thus 5(0) is given by 15(0) = Jgs Tr(a(o)(ia(o)). If : — > M is such that d{T,i) = ji, i.e. Sj 
a Seifert surface for ji, then we have that a(o) = ^ Y17=i and we get 



jj=l •^'S'^ i,j=l 



5(0) is a linear combination of linking numbers, so it is a link invariant as predicted from our 
general theory. We proceed to compute explicitly S'^i) which is given by 



5(1) = ^^Tr(2a(o)da(i) + 
We know that a(o) = ^ Y17=i and 



(o)J- 



da 



1 1 " 

-^Ko)>a(o)] + -E-^(7i'^i.iW)^' 



(1) 2'-"W"WJ ■ 2 



where ci = ^ X^^Li /ai jP(Ej)ci. Plugging these identities in the previous expression for 
5(1) we obtain 



i,j,k \ -^^^ -^^0,1 J 

The first summand in the formula above should be clear. The second summand arises from 

^Y.[ rr(P(S,)c,P(7„ et,,.(P(S,)cfe)c,), 

or equivalently 

By the defining properties of Poincare forms and antisymmetry of the Lie bracket, the later 
expression is equal to 

i,j,k 

The formula obtained for 5(i) is a link invariant with a crystal clear geometric meaning: 
the first summand counts triple intersections of the corresponding Seifert surfaces, the second 
summand counts pairs of points s, t in the parametrization of loop 7^, such that 7i(s) € and 
7i(t) € Tij. In the computation of 5(i) we make use of the identity 

n 

da{i) = -[a(o),a(o)] + ^ P{li,ci,i{t))ci, 



i=l 
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thus we assumed that the right hand side of this identity is a closed two-form. This assumption 
is by no means trivial and does not hold universally. Indeed it imposes a severe restriction 
on the type of links for which the invariant ^(i) is well-defined: the linking number of each 
pair of loops in the link must vanish. The Borromean rings is an example of link for which 
the invariant is well-defined and non-vanishing. For a proof of this and other interesting 
facts regarding the invariant 5(1) the reader may consult [23j. The reader should notice that 
^(i) is the second Milnor's invariant [25] for links in S^, and thus our method provides an 
interpretation for that invariant coming from perturbative Lagrangian physics. We expect that 
the higher order invariants S'(„) correspond with higher order Milnor's invariants which can be 
computed using higher order Massey products [26j . 



6 Yang-Mills theory and area invariants 

In this section we show that it is possible to obtain invariants of configurations of immersed 
curves in the plane from Yang-Mills- Wong action. To our knowledge results of this type have 
seldom been reported - unlike the relation between links and Chern-Simons theory - perhaps 
because the space of immersed curves in the plane, considered up to area preserving diffeomor- 
phisms, has not been deeply studied in the mathematical literature. The example consider in 
this section is studied in full details in [13] , here we only highlight the results of that paper that 
are useful to illustrate yet another application of our method. 

The basic settings is quite similar to those for Chern-Simons- Wong action. The space of 
fields is (g) s) x M(S'^,G)". The space of parameters /„(S'^,M^) consists of n-tuples 

(71, 7ra) such that 7^ : — > is an immersed closed curve in M?, and the images of the 7^ 
intersect, if they do, in transversal double points. The group of symmetries for the Yang-Mills- 
Wong action is the group of area preserving diffeomorphisms of M^. As before we fix q G g and 
for 1 < i < n we let Qi G M{S^, G). The action functional is given by 

S(a,5i,...,ffn,7i,-,7n) = [ Tr{FaA*Fa) + y2 [ dTTr{k,g7\T)Drgi{T)), 

where Fa = da + ^[a, a] and * is the Hodge star operator. According to our general theory the 
on-shell action may be expanded as 

00 

Sos{ll,-,ln) = ^ 'S'(„)(7i,...,7„)A" 
m=0 

where each S'^^) should be a function of (71, 7„) € /„(S'^, M^) invariant under area preserving 
diffeomorphisms of M^. One can compute S^q) ^i^d /S'(i) and check that they are indeed invariants 
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under area preserving diffeomorphisms. The invariant 5(o) admits the fairly simple expression 

where the functions J{"fi,^j) have the following geometric interpretation. A generic immersed 
curve in induces a partition of into a finite number of compact blocks and an unbounded 
block. The function J(7i,7j) is the sum of the signed areas of the intersections of the finite 
blocks of 7j with the finite blocks of 7j. In complete analogy with the Chern-Simons-Wong 
case, the geometric interpretation of 5(1) takes into account no just the areas of the intersection 
blocks, but also the order in which the intersection blocks appear for several, at least three, 
curves. Again 5(1) is only well-defined for an appropriated choice of curves. In p3] we describe 
explicitly three planar curves - a planar version of the Borromean rings - for which S'^j^) is 
well-defined and non- vanishing. 

7 Final remarks 

We introduced a method that yields invariant functions from classical field theories. In the 
perturbative regime we actually obtain a countable hierarchy of invariants. Our construction 
leaves many open problems and suggest new lines of research. For Chern-Simons-Wong action 
and 2-dimensional Yang-Mills- Wong action we are, at this point, only able to compute the first 
two invariants of the hierarchy. Though Theorem [16] provides explicit formulae for the higher 
order invariants, and our computations suggest that the consistency equations are satisfied, a 
rigorous proof is needed. We expect the higher order link invariants arising in the computation 
of the perturbative Chern-Simons-Wong on-shell action, to be closely related to Milnor's link 
invariants |25|. 

We believe our methods can be usefully applied to other classical field theories. In partic- 
ular, it may be rewarding to look at Yang-Mills- Wong action in higher dimensions, it should 
yield conformal invariants associated with closed curves in spacetime. It may also be interesting 
to apply our methods to the generalized Chern-Simons action of [29], it should yield invariants 
related with Chas-Sullivan product in string topology [9], [TO]. In our study of perturbative solu- 
tions we saw that the invertibility of the quadratic part of the action plays a fundamental role. 
Often the quadratic part is not invertible and new techniques are required in order to get invari- 
ants. One possibility is to introduce, as in the quantum case, fermionic variables and replace the 
action with a new one with invertible quadratic part. Thus, it is plausible that in the classical 
perturbative regime, the BRST and BV procedures may still play a role. Another possibility 
arises when the inverse of the quadratic part of the action is no quite well-defined, but rather a 
singular operator. In this case techniques from renormalization [H] may become useful in order 
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to replace invariants given by ill-defined divergent integrals, by their renormalized values. In 
recent years it has become clear that many constructions in field theory [2|I121[T^. as well as in 
other branches of physics and mathematics O El El [3 [2T] admit categorical analogues. It would 
be interesting to investigate the categorical foundations of the method introduced in this paper. 
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